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Abstract We present a general framework to calcu- 
late the properties of relativistic compound systems 
from the knowledge of an elementary Hamiltonian. Our 
framework provides a well-controlled nonperturbative 
calculational scheme which can be systematically im- 
proved. The state vector of a physical system is calcu- 
lated in light-front dynamics. From the general proper- 
ties of this form of dynamics, the state vector can be 
further decomposed in well-defined Fock components. 
In order to control the convergence of this expansion, we 
advocate the use of the covariant formulation of light- 
front dynamics. In this formulation, the state vector 
is projected on an arbitrary light-front plane uj-x = 
defined by a light-like four- vector uj. This enables us 
to control any violation of rotational invariance due to 
the truncation of the Fock expansion. We then present 
a general nonperturbative renormalization scheme in 
order to avoid field-theoretical divergences which may 
remain uncancelled due to this truncation. This general 
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framework has been applied to a large variety of mod- 
els. As a starting point, we consider QED for the two- 
body Fock space truncation and calculate the anoma- 
lous magnetic moment of the electron. We show that it 
coincides, in this approximation, with the well-known 
Schwinger term. Then we investigate the properties of 
a purely scalar system in the three-body approxima- 
tion, where we highlight the role of antiparticle degrees 
of freedom. As a non-trivial example of our framework, 
we calculate the structure of a physical fermion in the 
Yukawa model, for the three-body Fock space trunca- 
tion (but still without antifermion contributions). We 
finally show why our approach is also well-suited to 
describe effective field theories like chiral perturbation 
theory in the baryonic sector. 
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1 Light-front dynamics in few-body systems 
and field theory 

The understanding of hadron properties from an un- 
derlying Lagrangian or Hamiltonian is a major issue in 
nuclear and particle physics. It demands both a rela- 
tivistic framework to deal with quasi-massless particles 
(the pion, up and down quarks, etc.) or with high mo- 
mentum and high energy experiments, and a nonper- 
turbative framework. The latter is mandatory in order 
to calculate for instance the mass of a bound state from 
the pole of the scattering amplitude or from an eigen- 
state equation. 

One may already gain some physical insights from 
nonrelativistic studies like the nonrelativistic constituent 
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quark model in particle physics, or the study of few- 
nucleon systems, based on the nonrelativistic nucleon- 
nucleon or three-nucleon potentials, in nuclear physics. 
This has given rise to numerous studies in the last 40 
years. 

The extension to relativistic calculations may either 
rely on the use of relativistic equations like the four- 
dimensional Bethe-Salpeter equation (including its var- 
ious three-dimensional quasipotential reductions) and 
the Dyson- Schwinger equation, or on Hamiltonian dy- 
namics, using one of its three forms proposed by Dirac 
in 1949 pp. We shall follow in this review the path 
pointed out by Dirac and choose light-front dynamics 
(LFD) as a basis of our approach. 



1.1 Few-body relativistic systems 

A natural testing ground for the use of LFD is the study 
of few-body systems. The properties of this particular 
form of dynamics are indeed very well-suited to make 
a tight connection with nonrelativistic considerations. 

In the standard form of LFD, the state vector of 
a physical system is defined not at a fixed moment of 
time but on the light-front plane given by the equa- 
tion t + | = const. The nonrelativistic limit reached 
by taking c — > oo leads thus naturally to the ordinary 
equal-time formulation t = const, giving rise, in partic- 
ular, to the Schrodinger equation for the nonrelativistic 
wave function. 

One can always decompose the state vector in Fock 
components. Since the physical vacuum is trivial in 
LFD, i.e. it coincides with the vacuum for not interact- 
ing particles, this decomposition does not include the 
vacuum fluctuations but contains the physical states 
only. Few-body systems are represented just by the first 
few components of this expansion. Moreover, each Fock 
component has a probabilistic interpretation similar to 
that of the nonrelativistic wave function. 

In this formulation, the elementary kernel of the 
eigenstate equation defining the state vector of the phys- 
ical system is calculated in light-front time-ordered per- 
turbation theory. All four-momenta are on the mass 
shell, while all intermediate states are off the energy 
shell. The eigenstate equation to be solved is then three- 
dimensional, in direct analogy with the nonrelativistic 
Schrodinger equation. 

According to the Dirac's classification, the ten gen- 
erators of the Poincare group, given by space-time trans- 
lations (four generators), space rotations (three gener- 
ators), and Lorentz boosts (three generators), can be 
separated into kinematical and dynamical operators. 
The kinematical operators leave the light-front plane 



invariant and are independent of dynamics, i.e. of the 
interaction Hamiltonian of the system, while the dy- 
namical ones change the light-front position and depend 
therefore on the interaction. Among the kinematical op- 
erators, one finds, in LFD, the boost along the z axis. 
This property is of particular interest when one cal- 
culates electromagnetic observables at high momentum 
transfer, since once one knows the state vector in one 
reference frame, it is easy to calculate it in any other 
frame. 

One has of course to pay some price for that. The 
spatial rotations in the xz and yz planes become dy- 
namical, in contrast to the case of equal-time dynam- 
ics. This is a direct consequence of the violation of ro- 
tational invariance caused by the non-invariant defini- 
tion of the light-front plane orientation. This violation 
should be kept under control. 

Control of rotational invariance. While rotational in- 
variance should be recovered automatically in any ex- 
act calculation, this is not a priori the case if the Fock 
expansion is truncated. The control of the violation of 
rotational symmetry is very difficult in practice, when 
using the standard form of LFD. To avoid such an un- 
pleasant feature of the latter, we shall use below the 
covariant formulation of LFD (CLFD) [2j|3], which pro- 
vides a simple, practical, and very powerful tool in order 
to describe physical systems as well as their electromag- 
netic amplitudes. In this formulation, the state vector 
is defined on the plane characterized by the invariant 
equation uj-x = 0, where uj is an arbitrary light-like 
(uj 2 = 0) four-vector. The standard LFD on the plane 
t + | = is recovered by considering the particular 
choice uj — (1,0, 0,-1). The covariance of our approach 
relies on the invariance of the light-front plane equation 
under any Lorentz transformation of both uj and x. This 
implies in particular that uj cannot be kept the same in 
any reference frame, as it is the case in the standard 
formulation of LFD. 

There is of course equivalence, in principle, between 
the standard and covariant forms of LFD in any ex- 
act calculation. Calculated physical observables must 
coincide in both approaches, though their derivation in 
CLFD in most cases is much simpler and more trans- 
parent. Indeed, the relation between CLFD and stan- 
dard LFD reminds that between the Feynman graph 
technique and old-fashioned perturbation theory. 

In approximate calculations however, CLFD has a 
definite advantage in the sense that it enables a direct 
handle on the contributions which violate rotational in- 
variance. These ones depend explicitly on the orienta- 
tion of the light-front surface (i.e. on uj) and can thus be 
separated covariantly from true physical contributions. 
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This is of particular interest when one considers for 
instance electromagnetic observables in few-body sys- 
tems [3 j, or field theory on the light front. 

1.2 Light-front field theory 

The interest to the application of LFD to field-theoretical 
problems is also very old. It originates from the study 
of deep inelastic scattering experiments. Indeed, it is 
easy to see that in this kinematical domain, all events 
are close to the light-front plane. This legitimates the 
use of the infinite momentum frame in the calculation 
of high energy observables [4. . The calculation is done 
in perturbation theory. The extension to QED is natu- 
ral, since perturbation theory is also applicable in that 
case [5]. 

The use of LFD in field theory is important in or- 
der to extend these calculations to the nonperturbative 
domain [6]. As already mentioned, the structure of the 
vacuum in LFD enables a well-defined expansion of the 
state vector in Fock components. In field theory, this ex- 
pansion is in principle infinite. From a practical point 
of view, it should be truncated to a limited number of 
components. We thus have to worry about the conver- 
gence of this expansion, and the various ways to speed 
it up if one wants to be able to make meaningful pre- 
dictions. 

The truncation of the Fock expansion induces how- 
ever two pernicious features in the study of field theory 
on the light front. The first one, which we have already 
addressed in the discussion of few-body systems, is the 
violation of rotational invariance due to the particular 
choice of the orientation of the light-front plane. The 
second one is the appearance of uncancelled divergences 
which calls for an appropriate renormalization scheme. 

Appropriate renormalization scheme. The truncation of 
the Fock expansion complicates the renormalization pro- 
cedure, in contrast to that in standard perturbation 
theory. Indeed, the full cancellation of field-theoretical 
divergences which appear in a given Fock sector requires 
taking into account contributions from other sectors. If 
even a part of the latter is beyond our approximation, 
some divergences may leave uncancelled. 

For instance, looking at Fig. [T] for the calculation of 
the fermion propagator in the second order of pertur- 
bation theory, one immediately realizes that the cancel- 
lation of divergences between the self-energy contribu- 
tion (of 2nd order in the Fock decomposition) and the 
fermion mass counterterm (of 1st order one) involves 
two different Fock sectors. 

This means that, as a necessary condition for the 
cancellation of divergences, any mass counterterm should 
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Fig. 1 Renormalization of the fermion propagator in the second 
order of perturbation theory. The fermion mass counterterm is 
denoted by 5m. 



be associated with the number of particles present (or 
"in flight") in a given Fock sector. In other words, all 
mass counterterms must depend on the Fock sector un- 
der consideration, as advocated first in Ref. [7]. This is 
also true for the renormalization of the bare coupling 
constant. 

The presence of uncancelled divergences reflects it- 
self in possible dependence of approximately calculated 
observables on the regularization parameters (e. g., cut- 
offs). In other words, calculated physical observables are 
not anymore scale invariant. This prevents to make any 
physical predictions if we cannot control the renormal- 
ization procedure in one way or another. We have devel- 
oped an appropriate renormalization procedure — the 
so-called Fock sector dependent renormalization (FSDR) 
scheme — in order to keep the cancellation of field- 
theoretical divergences under permanent control. This 
scheme relies also directly on CLFD in order to define 
the renormalization conditions imposed on true physi- 
cal observables [8]. 

We shall detail in the following the main features of 
CLFD, as well as FSDR, and apply our general frame- 
work to various physical systems. 

2 Covariant formulation of light-front dynamics 

The state vector, ^ a (p), of a compound system cor- 
responds to definite values for its mass M, its four- 
momentum p, and its total angular momentum J with 
the projection a onto the z axis in the rest frame, i.e., 
the state vector forms a representation of the Poincare 
group. It satisfies the following eigenstate equations: 

p P <Pi a (p)=p P <i>t a (p), (i) 

P 2 <p J J{p) = M 2 «^», (2) 

5 2 4> J J{p) = -M 2 J(J+l)<f> J J(p), (3) 

5 3 <t> J J(p) = M a<t> J J(p), (4) 

where P p is the four-momentum operator, S p is the 
Pauli-Lubanski vector 
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and J is the four-dimensional angular momentum op- 
erator which is represented as a sum of the free and 
interaction parts: 

Jpv — JpJ + J pv • (6) 
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In terms of the interaction Hamiltonian H int {x) we 
have 

Jfv = / H int (x)(x p uj u - x u uj p )S(uj-x) d 4 x. (7) 

Similarly to J, the momentum operator also can be 
split into the free and interaction parts: 



x p ~ 
with 



p(0) 

p 



uj p / H int (x) S(u-x)d 4 x. 



(8) 



(9) 



From the general transformation properties of both 
the state vector and the light-front plane, it follows [9] 
that 

>(p) = VM^( p ), (io) 



where 



L pv (uj) 



d 



d 



(11) 



The equation (10) is called the angular condition. We 
can use it in order to replace the operator J™* entering 
into Eq. ^ by L pu (uj). Introducing the notations 



pv 



Wo 



7(0) 
d pv 

1 



, £pva/3 



P v M 



aft 



(12) 
(13) 



we obtain, instead of Eqs. (|3| and Q, 



W 2 4> j J{p) = -MV(J+i) <p J J(p), 



Jal 



(14) 

W 3 4>i a (p) =Ma 4>i a (p)- (15) 
These equations do not contain the interaction Hamil- 
tonian, once (j)^f satisfies Eqs. ([!]) and The con- 
struction of the state vector of a physical system with 
definite total angular momentum becomes therefore a 
purely kinematical problem. Indeed, the transformation 
properties of the state vector under rotations of the co- 
ordinate system are fully determined by its total angu- 
lar momentum, while the dynamical part of the latter is 
separated out by means of the angular condition. The 
dynamical dependence of the state vector on the light- 
front plane orientation turns now into its explicit depen- 
dence on the four- vector uj [3]. Such a separation, in a 
covariant way, of kinematical and dynamical transfor- 
mations is a definite advantage of CLFD, as compared 
to standard LFD. 



2.1 General Fock decomposition of the state vector 

According to the general properties of LFD, we decom- 
pose the state vector of a physical system in Fock sec- 
tors. We have 

OO r. 

= / dD n(t>n(k 1 ,...,k n ]p) 



x £ 4 (/ci + • • • + K - p - ur n ) \n) , 



(16) 



where |n) is the state containing n free particles with 
the four- momenta fci, . . . , k n and n 's are relativistic n- 
body wave functions or the so-called Fock components. 
Here and below we will omit, for shortness, all spin in- 
dices in the notation of the state vector. Note the partic- 
ular overall momentum conservation law given by the 
J-function. It follows from the general transformation 
properties of the light-front plane uj-x = under four- 
dimensional translations. The quantity r n is a measure 
of how far the n-body system is off the energy shelf] It 
is completely determined by this conservation law and 
the on-mass-shell condition for each individual particle 
momentum. We get 



2uj-pr n = (s n -M 2 ), 
where 

s n = (h + . . . + k n ) 2 . 



(17) 



(18) 



The phase space volume element is represented schemat- 
ically by dD n . It involves integrations over the compo- 
nents of the constituent four- moment a d 4 ki and over 
dr n in infinite limits. The state |n) can be written as 

|n) =d t (fci)d t (fe)...d t (fc„)|0) ) (19) 

where d^ is a generic notation for free particle creation 
operators. To completely determine the state vector, we 
normalize it according to 



0(p') t 0(p) = 2po^ 3) (p'-p)- 



(20) 



It is convenient to introduce, instead of the wave func- 
tions n , the vertex functions r n (which we will also 
refer to as Fock components), defined by 



(s n - M 2 )(j) n 2uj-p r n <j) n 



(21) 



In the particular case of a fermion coupled to bosons, it 
is convenient to extract from r n the fermion bispinors, 
and make the replacement 



(22) 



where k\ is the four-momentum of the constituent fermion. 
When Fock space is truncated, it is necessary to keep 
track of the order of truncation N (i.e. the maximal 
number of particles admitted in the Fock sectors) in 
the calculation of the vertex function. For this pur- 
pose we will use the notation for the n-body ver- 
tex function. In the LFD graph technique, it is repre- 
sented by a (n-hl)-leg vertex with one incoming double 
line corresponding to the physical state and n outgoing 
single lines corresponding to constituents. By its spin 
structure and transformation properties it is completely 

1 The term "off the energy shell" is borrowed from the 
equal-time dynamics where the spatial components of the four- 
momenta are always conserved, but the energies of intermediate 
states are not equal to the incoming energy. 
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analogous to a n-body wave function. As an example, 
we show such a vertex in Fig. [2] for the case of a physical 
fermion state composed from one constituent fermion 
and (n — 1) bosons. For simplicity, we shall use be- 
low the same notation for the vertex functions of both 
fermion and boson physical states. 



(n — 1) bosons 




Fig. 2 n-body vertex function for a Fock space truncation of 
order N, for a physical fermion (double straight line) made of a 
constituent fermion (single straight line) coupled to bosons (wavy 
lines). 



With the decomposition ( 16 ), the normalization con- 



dition (20) writes 

oo 

£/„ = !, 



(23) 



where I n is the contribution of the n-body Fock sector 
to the full norm of the state vector, equal to 1. The 
explicit formulas for I n in terms of the vertex functions 
for some important particular cases can be found in [8] . 



2.2 Eigenstate equation 

The system of coupled equations for the Fock compo- 
nents of the state vector can be obtained from Eq. ([2| 
by substituting there the Fock decomposition ( fl6| ) and 
calculating the matrix elements of the operator P 2 in 
Fock space. With the expressions ([8| and ([9|, we get 
the eigenstate equation [TP] : 

dr 



l(u>-p) J 



(pv»y 



M 



Hp), 

(24) 

where H mt is the interaction Hamiltonian in momen- 
tum space: 



H in \oj T ) = / H int (x)e-^' x)T d 4 x. 



(25) 



Using the general momentum conservation law in 



Eq. (16), we conclude that the operator in the square 



brackets on the right-hand side of Eq. ( 24 ) simply mul 



tiplies each Fock component of the state vector by the 
factor s n — M 2 = 2(u-p)r n . It is therefore convenient to 
introduce the notation 



where f is the operator which, acting on a given compo- 
nent (j) n of 0(p), gives r n (j) n . Q(p) has a Fock decomposi- 



tion which is obtained from Eq. ( 16 ) by the replacement 
of the wave functions (j) n by the vertex functions r n . We 
can thus cast the eigenstate equation in the form 



G(p) 



l 

2tt 



dr 

-H mt (ujT) —G(p) 



(27) 



The physical mass M of the compound system is found 
from the condition that the eigenvalue is 1. This equa- 
tion is quite general and equivalent to the eigenstate 
equation It is nonperturbative. 



2.3 Spin decomposition of the state vector 

As follows from the angular condition, the spin struc- 
ture of the wave functions (j) n is very simple, since its 
construction does not require the knowledge of dynam- 
ics. It should incorporate however cj-dependent compo- 
nents. It is convenient to decompose each wave function 
(j) n into invariant amplitudes constructed from all avail- 
able particle four-momenta (including the four-vector 
col) and spin structures (matrices, bispinors, etc.). In 
the Yukawa model for instance, we have for the one- 
and two-body components [8]: 



01 = ipi u(ki)u(p), 



<j>2 u(ki) 



uj-p 



u{p), 



(28) 
(29) 



g(p) = 2(u-p)f(j)(p), 



since no other independent spin structures can be con- 
structed. Here ^i, ^2, and ^ 2 are scalar functions de- 
termined by dynamics. For a spin 1/2 physical fermion 
composed from a constituent spin 1/2 fermion coupled 
to scalar bosons, the number of invariant amplitudes for 
the two-body Fock component coincides with the num- 
ber of independent amplitudes of the reaction spin 1/2+ 
scalar — >• spin 1/2 + scalar, which is (2 x 2)/2 = 2, due 
to parity conservation. Similar expansions can be done 
for the three-body component of the same system or for 
Fock components in QED, as we shall see in Sec. [4j 



3 Fock sector dependent renormalization 
scheme 

In the standard renormalization theory, the bare pa- 
rameter^] are determined by fixing some physical quan- 
tities like the particle masses and the physical coupling 
constants. The bare parameters are thus expressed thro- 
ugh the physical ones. This identification implies in fact 

2 The term "bare parameters" means here the whole set of 
parameters entering into the interaction Hamiltonian, e.g. the 



(26) bare coupling constants and the mass counterterms. 
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the following two important questions which are usually 
never clarified in LFD calculations, but are at the heart 
of our scheme. 

(i) In order to express the bare parameters through 
the physical ones, and vice- versa, one should be able to 
calculate observables or, in other words, physical ampli- 
tudes. In LFD, any physical amplitude is represented as 
a sum of partial contributions, each depending on the 
light-front plane orientation. Since an observable quan- 
tity can not depend on the latter, this spurious depen- 
dence must cancel in the whole sum, as already men- 
tioned. Such a situation indeed takes place, for instance, 
in perturbation theory, provided the regularization of 
divergencies in LFD amplitudes is done in a rotation- 
ally invariant way [TT]. In nonperturbative LFD calcu- 
lations, which are always approximate, the dependence 
on the light-front plane orientation may survive even 
in calculated physical amplitudes. For this reason, the 
identification of such amplitudes with observable quan- 
tities becomes ambiguous and expressing the calculated 
amplitudes through the physical parameters turns into 
a non-trivial problem. 

(ii) The explicit form of the relationship between 
the bare and physical parameters depends on the ap- 
proximation which is made. This is a trivial statement 
in perturbation theory where the order of approxima- 
tion is distinctly determined by the power of the cou- 
pling constant. In our nonperturbative approach based 
on the truncated Fock decomposition of the state vec- 
tor, an analogous parameter (the power of the coupling 
constant) is absent. At the same time, to make cal- 
culations compatible with the order of truncation, one 
has to trace somehow the level of approximation. This 
implies that, on general grounds, the bare parameters 
should depend on the Fock sector in which they are 
considered. Moreover, this dependence must be such 
that all divergent contributions are cancelled, as already 
mentioned in Sec. 11.21 How this should be done is the 
crucial point of the nonperturbative formulation of field 
theory on the light-front. 

The use of our FSDR scheme in CLFD is a unique 
opportunity to answer both questions. For clarity, we 
shall take, as a background, a model of interacting fer- 
mions and bosons like the Yukawa model or QED, with 
the aim to calculate the physical fermion state vector. 
The basis of Fock space is formed by a set of Fock sec- 
tors, each containing one constituent fermion and a cer- 
tain number of bosons (0, 1, 2,...). The truncation is 
made by retaining only those Fock sectors where the 
maximal total number of constituent particles does not 
exceed N. The consideration of antiparticle degrees of 
freedom will be discussed in the simple case of a scalar 
model in Sec. 14.31 We will assume that the interaction 



Hamiltonian H int {x) is constructed through the bare 
fields satisfying the free Dirac or Klein- Gordon equa- 
tions with the corresponding physical masses, while the 
mass renormalization is performed by introducing, into 
H tnt (x), appropriate mass counterterms. We emphasize 
that in spite of that we consider hereafter several par- 
ticular forms of interaction, our scheme is applicable to 
physical systems with arbitrary interaction admitting a 
Fock decomposition of the state vector. 

3.1 Mass counterterm 

The simple example of the renormalization of the fermion 
self-energy within the two-body Fock space truncation, 
presented in Sec. |1.2[ can serve as a guideline to set up 
our general rule. In this example, the mass countert- 
erm should be labelled with a subscript and denoted 
by Jm2, in order to indicate that it is introduced to 
cancel, at $ = m, where m is the constituent fermion 
mass, the self-energy contribution which belongs to the 
two-body Fock sector. Let us denote by Smi the mass 
counterterm in the most general case. Since we truncate 
our Fock space to order TV, one should make sure that, 
at any light-front time, the total number of particles is 
at most N. Our first rule is thus: 

— in any amplitude where the mass counterterm Smi 
appears, the value of I is such that the total number 
of bosons in flight plus I equals the maximal number 
of the Fock sectors considered in the calculation, i.e. 
N. 

For instance, in the typical contribution indicated in 
Fig. [31 the mass counterterm is Jm^-n+i)- For the 




(n — 1) bosons 



if Sm {N - n+1 ) 



Fig. 3 Typical insertion of the mass counterterm. 

mass counterterm of the lowest order, we simply have 

Sm 1 = 0, (30) 

since the fermion mass is not renormalized at all if the 
fermion can not fluctuate in more than one particle! 

3.2 Bare coupling constant 

The general strategy we developed above for the calcu- 
lation of the mass counterterm should be also applied 
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to the calculation of the bare coupling constant, with 
however a bit of caution, since this one may enter in 
two different types of contributions. 

The first one appears in the calculation of the state 



vector itself, when Eq. (27) is solved. In that case, any 



boson-fermion coupling constant is associated with the 
emission or the absorption of a boson which participates 
in the particle counting, in accordance with the rules 
detailed above, since it is a part of the state vector. 

The second one appears in the calculation of the 
boson-fermion scattering amplitude or of the boson- 
fermion three-point Green's function (3PGF) like the 
electromagnetic form factor. Since the external boson 
is an (asymptotic) free field rather than a part of the 
state vector, the particle counting rule advocated above 
should therefore not include the external boson line. 

Following the same reasoning developed above for 
the calculation of mass counterterms, we can formulate 
the following general rule for the calculation of the bare 
coupling constants 

— in any amplitude which couples constituents inside 
the state vector one should attach to each vertex the 
internal bare coupling constant goi • The value of I 
is such that the total number of bosons in flight be- 
fore (after) the vertex - if the latter corresponds to 
the boson emission (absorption) - plus I equals the 
maximal number of the Fock sectors considered in 
the calculation, it i.e. N. 

The calculation of external bare coupling constants pro- 
ceeds in the same spirit, with the final rule: 

— in any amplitude which couples constituents of the 
state vector with an external field, one should at- 
tach to the vertex involving this external field the 
external bare coupling constant g$i . The value of I 
is such that, at the light-front time corresponding to 
the vertex, the total number of internal bosons in 
flight (those emitted and absorbed by particles en- 
tering the state vector) plus I equals the maximal 
number of the Fock sectors considered in the calcu- 
lation, i.e. N. 

The lowest order bare coupling constants are 

9oi = , 0oi 9- (31) 

The first one is trivial, because no fermion-boson inter- 
action is allowed in the one-body Fock space truncation. 
The second one reflects the fact that the external bare 
coupling constant, in the same approximation, is not 
renormalized at all, since a single fermion can not be 
" dressed" . 

Some illustrations of the rules concerning the inter- 
nal and external bare coupling constants are given in 
Fig. [4} 



(a) 




| (n — 2) bosons 



90(N-n+2) 



(b) 




(n — 1) bosons 



Y { n ] 90(N-n+l) 



(c) 



gO(N-n+l) 



J- n 




-L n 



(n — 1) bosons 

Fig. 4 Typical contributions to the fermion state vector for the 
absorption (a) and the emission (b) of an internal boson, and to 
the fermion-boson 3PGF (c). 



Though we relied on the fermion-boson model when 
considering the above FSDR procedure, the latter can 
be easily extended to other systems with additional 
counterterms and bare parameters. 



3.3 Renormalization conditions and wave function 
renormalization 

Once proper bare coupling constants and mass coun- 
terterms have been identified, one should fix them from 
a set of renormalization conditions. In perturbation the- 
ory, there are three types of quantities to be deter- 
mined: the mass counterterms, the bare coupling con- 
stants, and the norms of the fermion and boson fields. 
Usually, the on-mass-shell renormalization is applied, 
with the following conditions. For each field, the mass 
counterterm is fixed from the requirement that the cor- 
responding two-point Green's function has a pole at 
p 2 = m 2 , where m is the physical mass of the parti- 
cle. The field normalization is fixed from the condition 
that the residue of the two-point Green's function at 
the pole is 1. The bare coupling constant is determined 
by requiring that the on-mass-shell 3PGF is given by 
the product of the physical coupling constant and the 
elementary vertex. 
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The renormalization conditions in LFD are of slightly 
different form, although they rely on the same grounds. 
The fermion mass counterterm is fixed from the eigen- 
value equation Q by demanding that the mass, M, of 
the physical state is identical to the mass of the con- 
stituent fermion which we called m. Solving a similar 
equation for the physical boson state vector allows us 
to find the boson mass counterterm from quite analo- 
gous requirements. The determination of the internal 
bare coupling constant needs some care. It is found, as 
in perturbation theory, by relating the on-energy-shell 
two-body vertex function F 2 to the physical coupling 
constant g. As follows from the momentum conserva- 
tion law, taking F 2 on the energy shell is equivalent to 
setting r 2 0, or s 2 = (fci + k 2 ) 2 = M 2 . 

In order to fix the relationship between F 2 and g one 
needs to take into account the renormalization factors 
coming from radiative corrections to all legs of the two- 
body vertex function [12] . These factors do also depend 
on the order of the Fock space truncation, as detailed 
in [15]. In the case of a fermion coupled to bosons, with- 
out polarization effects generated by antifermion con- 
tributions, this relationship reads 



r 2 w ( S2 = M 2 ) 



(N-l) 



(32) 



where I\ is the one-body contribution to the norm of 



the state vector, as given in Eq. (23), calculated for the 
Fock space truncation of order N — 1. 



Eq. (32) admits simple physical interpretation. Each 



leg of the on-energy-shell two-body vertex function con- 
tributes for an individual factor \[Z to the physical cou- 
pling constant, where Z is the field strength normaliza- 
tion factor [15] . The physical fermion state is normal- 
ized to 1, so that its factor Z equals 1. The constituent 
boson line is not renormalized - since we do not con- 
sider antifermions - so that for the bosonic line we also 
have Z\) — 1. If we did not neglect the antifermionic 
degrees of freedom, it would contribute by a non-unity 
factor \fZ\j. Finally, we have shown in [15] that the 
field strength normalization factor of the constituent 
fermion is just the weight of the one-body component 
in the norm of the physical state, i.e. Z$ = I\. Accord- 
ing to our FSDR scheme, the normalization factor of 
the constituent fermion should correspond to the Fock 
space truncation of order N — 1, since, by definition, 
the two-body vertex function contains one extra boson 
in flight in the final state. 



The condition (32) has two important consequences. 



The first one is that the two-body vertex function at 
s 2 = M 2 should be independent of the four-vector uj 
which determines the orientation of the light-front plane. 



With the spin decomposition (29), this implies that the 
component at s 2 = M 2 should be identically zero. 



While this property is automatically verified in the case 
of the two-body Fock space truncation - if using a regu- 
larization scheme which does not violate rotational in- 
variance - this is not guaranteed for calculations within 
higher order truncations. 

Indeed, nothing prevents F 2 to be cj-dependent, since 
it is an off-shell object, but this dependence must com- 
pletely disappear on the energy shell, i.e. for s 2 = M 2 . 
It would be so if no Fock space truncation occurs. The 
latter, in approximations higher than the two-body one 
(i. e. for N = 3, 4, . . .), may cause some cj-dependence 
of r 2 even on the energy shell, which immediately makes 



the general renormalization condition (32) ambiguous. 
If so, one has to insert new counterterms into the light- 
front interaction Hamiltonian, which explicitly depend 
on uj and cancel the cj-dependence of r 2 (s 2 — M 2 ). 
Note that the explicit covariance of CLFD allows to sep- 
arate the terms which depend on the light-front plane 
orientation (i.e. on uj) from other contributions and es- 
tablish the structure of these counterterms. This is not 
possible in ordinary LFD. 



We should thus enforce the condition ( 32 ) by intro- 



ducing, into the interaction Hamiltonian, appropriate u 
dependent counterterms. For instance, in the Yukawa 
model without antifermion contributions and within 



the three-body Fock space truncation (see below, Sec. 4.4) 
we need one additional counterterm of the form 



6HZ 



— Tfi (/) 

-z^—^^y, 

iuj-o 



(33) 



where Z^ is a constant adjusted to make Eq. (|32j) true, 
i/)((p) is the fermion (boson) field, and l/i(uj-d) is the 
reversal derivative operator, fully analogous to the op- 
erator l/id + in ordinary light-front dynamics. 



The second consequence of Eq. ( 32 ) is that the two- 



body vertex function at 52 = M 2 should be a constant. 
This is a non-trivial requirement, since i~2, as well as the 



components ^ 2 and ^ 2 in Eq. (29), do depend on two 
invariant kinematical variables which are usually chosen 
as the longitudinal constituent momentum fraction, x, 
and the square of its transverse momentum, k^. The 
on-mass-shell condition 



s 2 



M 



(34) 



where \i is the constituent boson mass, fixes one variable 
(say, k^) only, while the second variable (x) remains ar- 
bitrary. Hence, under the condition ( |34| ), if we assume 
k^ fixed, r 2 should be independent of x. Again, this 
property is verified in the two-body Fock space trun- 
cation, since in this approximation our equations are 
equivalent to perturbation theory of order g 2 . It is not 
guaranteed for higher order calculations. We shall come 
back to this point in Sec. [4] In practice, we shall fix 
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i~2(s2 = M 2 ) at some preset value x* and verify that 
the physical observables are not sensitive to the choice 
of x*. 

To summarize, we can thus list the normalization 
conditions in CLFD, for a calculation done in a Fock 
space truncation of order 7V, without considering fermion- 
antifermion polarization corrections: 

— The mass counterterm drriN is fixed by solving the 
eigenstate equation (27) in the limit M — >> m. 



The state vector is normalized according to the stan- 
dard condition (2^ . 

The internal bare coupling constant goiv is fixed from 
the condition that the uj -independent part of the two- 
body vertex function at S2 = m 2 and at a fixed value 
of x, denoted by x* , is given by the right-hand side 



ofEq. (32). 



— The external bare coupling constant #oiv is fixed from 
the condition that the uj -independent part of the on- 
energy-shell 3PGF is proportional to the elementary 
vertex, with the proportionality coefficient being the 
physical coupling constant. 

— The uj -dependent counterterms in the Hamiltonian 
are fixed by the conditions that the uj- dependent parts 
of the two-body vertex function and the 3PGF, both 
taken on the energy shell, turn to zero. 

— The values of all bare parameters and counterterms 
for I < N are determined from successive calcula- 
tions within the 1, 2, ... TV Fock space truncations. 

4 Applications 



The explicit solution of the eigenvalue equation (27) re 



quires to define first a regularization scheme, since loop 
contributions are a priori divergent. We shall use in 
the following applications the Pauli-Villars (PV) regu- 
larization scheme, first applied to LFD in [17]. It has the 
nice feature of being rotationally invariant [TT] , and it 
can be implemented rather easily in calculations within 
the two- and three-body Fock space truncations [8] . Be- 
sides that, in this regularization scheme, all contact in- 
teractions inherent to LFD are absent. It however ne- 
cessitates to extend Fock space in order to embrace 
PV fermions and PV bosons on equal grounds with the 
physical particles. 



4.1 Self-energy of a fermion in the Yukawa model 

Before discussing the calculation of the properties of 
compound systems, it is instructive to look at the struc- 
ture of the fermion self-energy in the simple Yukawa 
model in CLFD. The results are very similar to those 
for QED [IT]. 



Since our formalism is explicitly covariant, we can 
write down immediately the general structure of the 
self-energy of a fermion with the off-shell four-momentum 
k (k 2 m 2 ). It writes 



U(}t) = A- 



m 



uj-k 



(35) 



where a = {ftyb — yb^)/A{u-k). The coefficients in this 
expansion depend on k 2 only. They are given by 
1 



^(fc 2 ) 
B(k 2 ) = 

C(k 2 ) 
Ci(fc 2 ) 



rTr TO)], 



m 



-Tr 



uj-k 



(36) 
(37) 

(38) 
(39) 



Auj-k 
1 

4m 

Tr [E(fi)a} 

In the two-body approximation, when U(f6) is entirely 
given by the loop diagram shown in Fig. [I] the coeffi- 
cient C\ is identically zero. The coefficient C should also 
be zero since the two-point Green's function should be 
equivalent to the one calculated in the four-dimensional 
Feynman approach, and is therefore independent of uj 
provided one uses a rotationally invariant regulariza- 
tion [11 j. Note that the coefficient C is not a priori 
chiral invariant in the sense that if the mass of the 
constituent fermion goes to zero, C does not vanish, 
in contrast to A and S, as it should. Using the PV reg- 
ularization scheme, A and B depend logarithmically on 
the PV boson mass. Without regularization, the coef- 
ficient C diverges quadratically at high momenta. It is 
however identical to zero when the PV regularization 
scheme is used with one PV fermion and one PV boson 
only. 

4.2 QED in the two-body approximation 

This simple case provides a good starting point to un- 
derstand how our general framework should be applied 
in practice. The eigenstate equation one has to solve 
is shown graphically in Fig. [5] This equation is writ- 
ten in accordance with the prescriptions of our FSDR 
scheme. Note the use of the Fock sector dependent mass 
counterterms and the bare coupling constants. The two- 
body vertex function writes in that case 



u(*i)[r£ e £(fc2)]u(p), 



(40) 



where is the polarization vector of the photon. To 
simplify notations, we omit the superscript "(2)" at 
rJf. We shall use in the following the Feynman gauge. 
The number of independent invariant amplitudes for 
this vertex function coincides with that for the reaction 
1/2 + 0— 1/2 + 1. However, one should take into ac- 
count that in the Feynman gauge the vector boson wave 
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function has four independent components. So, the to- 
tal number of invariant amplitudes is (2 x 2 x 4)/2 = 8. 
We choose the following set of invariant amplitudes [18] : 

m 2 ubuo^ 



1 2 



&17 M + 






m ^7 M 


D 2 




uo-p 




uo-p 


p» 
+65 — 


u K 

+ ^6— " 




m (Jop^ 


m 


m 




uo-p 



(uo-p) 2 

m i/jki 

1 • 

uo-p 



(41) 



In the two-body approximation, and using the PV reg- 
ularization scheme, we find 



h = 2e 02 mV'i, 
b 2 -8 = 0, 



(42) 
(43) 



where ipi is defined in Eq. (28). These components re 



fer to the physical ones. The ones associated with PV 
bosons and/or fermions can be found easily [8]. In this 



.12) 



r (2) 
1 2 



(2) 



5m 2 




With all these quantities, one can easily calculate the 
electromagnetic form factors given by the diagrams shown 
in Fig. [6j Note that since the state vector is normalized 
to 1, we find for the external bare coupling constant: 



eoi = e 02 = e. 



(51) 



At zero momentum transfer, the Pauli form factor gives 
the anomalous magnetic moment of the electron. We 

get i 

^2(0) = (52) 

This expression exactly coincides with the well-known 
Schwinger term. One may wonder why our nonpertur- 
bative approach exactly recovers the perturbative re- 
sult [T9] , The reason is simple. For the two-body Fock 
space truncation we consider here, the irreducible con- 
tributions to the two-point Green's function or to the 
electromagnetic form factors are just identical to the 
corresponding ones in the second order of perturbation 
theory. And their re-summation to all orders of the cou- 
pling constant just defines the physical mass of the elec- 
tron, which is also the same, by construction, in both 
approaches. 



Fig. 5 System of equations for the vertex functions in the two- 
body approximation. 



approximation, the vertex function F 2 is a constant 
matrix proportional to 7 M . With these results, one can 
calculate the norms of the Fock sectors entering into 



the normalization condition (23): 



J x = 4m 2 ^ 2 , 

h = 4ra 2 ^i e% 2 J 2 , 



(44) 
(45) 

where the expression for J2 can be found in [18]. It 
depends logarithmically on the mass of the PV boson 
used to regularize the loop integral. The normalization 
condition thus fixes ipi m . 
1 



1 



-02 



J2 



(46) 



The renormalization condition ( 32 ) enables us to calcu- 
late eo2 as a function of the physical coupling constant 
denoted by e. This condition, for N = 2, writes simply 

2e 02 rm/>i7 /i , (47) 



1 2 



since = 1. With Eq. (|46j) we get 



-02 



e 2 J 2 



and 

h = l-e 2 J 2 , 
h = e 2 J 2 . 



(48) 

(49) 
(50) 




TY> e 02 TY> 

Fig. 6 Fermion-boson 3PGF in the two-body approximation. 



Our result ( 52 ) is also independent of the P V fermion 
and boson masses, when they tend to infinity. This is 
of course necessary in order to preserve the scale invari- 
ance of physical observables. Note that at large enough 
PV masses, the one-body part I\ of the normalization 
condition is negative, while the two-body part I 2 ex- 
ceeds 1. This implies also that e\ 2 is negative. These 
features are just an artefact of the regularization scheme 
which is used. They show that indeed I\ and I 2l and, 



more generally, any I n in Eq. (23), as well as the bare 
coupling constants, are not physical observables and 
can therefore be scale dependent. 



4.3 Scalar system in the three-body approximation: 
the role of antiparticles 

Finding the state vector for the three-body Fock space 
truncation is the first non-trivial nonperturbative cal- 
culation. We begin our consideration of the three-body 
approximation with the study of a toy model: a heavy 
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scalar boson B with mass m, interacting with light 
scalar bosons b with mass fi. We will calculate the state 
vector of the heavy boson and represent the former as 
a sum of the following Fock sectors, in symbolic form: 

<t>{p) = \B) + \Bb) + \Bbb) + \BBB), (53) 

where B denotes the antiparticle. The interaction Hamil- 
tonian is 

H in \x) = -g B 2 (x)b(x) - Sm 2 B 2 (x) - £/iV(x), (54) 

where Sm 2 and 5ji 2 are the mass counterterms which, in 
this model, have the dimension of mass squared. Since 
the physical coupling constant g has the dimension of 
mass, it is convenient to define a dimensionless coupling 
constant a = g 2 /167rm 2 . 

The system of eigenstate equations for the vertex 
functions is represented graphically in Fig. [7J The two 
three-body vertices Jg 3 ^ and rj^ are trivially expressed 

(3) 

through the two-body vertex . Due to this fact, we 
can obtain a closed equation involving r£ only. 



p(3) 






r-(3) 
1 3 



1 3a 




Fig. 7 System of equations for the vertex functions of the scalar 
system in the three-body approximation. The thin and thick or- 
dinary lines stand for heavy bosons and antibosons, respectively. 
The wavy lines denote light bosons. 



In this model, we encounter two types of irreducible 
divergent contributions: the heavy and light boson self- 
energies, as shown in Fig. [8] Both of them diverge log- 
arithmically at high momenta; the divergences are ex- 
pected to be cancelled by the corresponding mass coun- 
terterms. 



Fig. 8 Divergent contributions in the scalar model: the heavy 
(a) and light (b) boson self energies. 



The renormalization condition should now incorpo- 
rate the light boson field strength normalization factor. 



Instead of Eq. ([32]) , we have 

ZBb, 



(3)/ 



2 {s 2 = M z )=g^Z Bh , (55) 

where the factor Zsb stands for the combined contri- 
bution of the heavy and light boson strength normal- 
ization factors, calculated according to the FSDR pre- 
scription, extended to include antiparticle degrees of 
freedom [13 . This factor is completely finite in a pure 
scalar system, and is expressed in terms of derivatives 
of the heavy and light boson self-energies. 

(3) 

The eigenstate equation for with the renormal- 



ization condition (55) has been solved numerically in 
the limit M — >• m for the following set of parameters: 
m = 0.95, a = 0.15, a = 2. Note that this value of a is 

(3) 

far from the perturbative domain. For convenience, r% 
is parameterized as a function of two variables, S2 and 
x. The physical domain corresponds to S2 > (m + /i) 2 , 
< x < 1. The renormalization point is defined by 
S2 = m 2 and x = x*, where we chose = /i/(ra + /i). 
We show in Fig. [9] the characteristic dependence of 

( 3^ 

i~2 on one °f ^s arguments, while the second argument 

(2) 

is fixed. In contrast to the two-body case, where r 2 ^ s 
a constant, now it exhibits rather nontrivial behavior, 
especially as a function of 52 at fixed x. 

To check the consistency ofthe renormalization con- 



dition (55), we show in Fig. 
on x at fixed 52 = m 2 
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the dependence of R 



It is distinctly seen that if the 



Fock space includes all the four Fock sectors, Eq. (53), 



rn 



,x) is almost a constant, as it ought to 
be. If one however neglects one of the three-body sec- 
tors, either \BBB) or \Bbb), = m 2 ,x) strongly 
depends on x, which makes the renormalization condi- 
tion ambiguous because physical results turn out to be 
sensitive to the choice of x*. 

The above result is very non-trivial. Let us consider, 
for instance, a part of the three-body contributions to 

(3) 

i~2 calculated in perturbation theory, i.e. in terms of 



the two-body vertex function R 



(2) 



two-body, approximation, as shown in Fig. 11 



found in the lower, 
Their 



sum (but not each of them!) is exactly ^-independent 
at S2 = m 2 , since is a constant. This statement 
follows from the fact that the sum of the two contribu- 
tions on the energy shell coincides with the correspond- 
ing on-mass-shell Feynman amplitude, and the latter is 
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an arbitrary chosen function i~2(s2jz0 instead of , 
their sum would hardly be a constant at S2 = m 2 . Such 
a result may indicate that the set of contributions we 
considered in the case of the three-body Fock space 
truncation forms an almost consistent set, in the sense 
that it meets all consistency requirements. This consti- 
tutes a basis of a well-controlled, nonperturbative Fock 
expansion of the the state vector, in the same spirit as 
an expansion in g 2 in perturbation theory. 




#02 



Fig. 11 Part of the three-body Fock sector contributions to the 
two-body vertex. 
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Fig. 9 Two-body vertex function in the scalar model as a func- 
tion of x at fixed S2 (upper plot) and as a function of S2 at fixed 
x (lower plot). 
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Fig. 10 Two-body vertex function (in units of gV^Bb) a t S2 = 
m 2 as a function of x. The solid line is the calculation with all the 
four Fock sectors, Eq. ( |53| ), in the state vector. The dash-dotted 
and dotted lines are the same but without the \BBB) or \Bbb) 
sector, respectively. 



a constant. In our calculation of order N — 3, the con- 
tributions analogous to those in Fig. 111! but with F^ 

(2) 

instead of appear in the eigenstate equation and 
determine the ^-dependence of /^(^ = rn 2 ). From 
Fig. [lO] it is seen that the latter dependence is surpris- 
ingly weak, though is not a constant, as Fig. [9] 
clearly demonstrates. At the same time, if one calcu- 



4.4 Yukawa model in the three-body approximation: 
calculation of the fermion anomalous magnetic 
moment 

The Yukawa model is much closer to reality (e.g. for 
instance to QED) than the scalar model discussed in 
the previous section. Simultaneously, it is much more 
involved from the point of view of renormalization, be- 
cause of complicated spin structure of the vertex func- 
tions and stronger divergences to be regularized and 
renormalized. In this section, we apply our FSDR renor- 
malization scheme to the Yukawa model within the three- 
body Fock space truncation, but, in order to avoid extra 
complications, without incorporating antifermion de- 
grees of freedom. Our purpose is to demonstrate the ca- 
pabilities of the FSDR scheme in solving a true nonper- 
turbative problem for particles with spin. For example, 
in contrast to the two-body Fock space truncation (in 
this approximation the Yukawa model is almost equiva- 



lent to QED, as considered in Sec. 4.2), where we have 
only one irreducible contribution to the fermion self- 
energy, we now should consider all the graphs for the 
self-energy, which contain one fermion and two bosons 
in intermediate states, including overlapping self-energy 
type diagrams. The number of such irreducible graphs 



is infinite. Some of them are shown in Fig. 12 The so- 
lution of the eigenstate equation for the state vector 
automatically generates these contributions to all or- 
ders of the coupling constant. 

We start with the following Fock decomposition for 
the state vector 



lates the amplitudes of the diagrams in Fig. 11 with cf)(p) = \F) + \Fb) + \Fbb), 



(56) 
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Fig. 12 Radiative corrections to the self-energy. 

where the symbols F and b denote, respectively, the 
constituent fermion and boson with masses m and ja. 
The interaction Hamiltonian reads 



(57) 



H in \x) = -g W<p - SrmjjiP - S^ 2 ip 2 + 5H™\ 
where SH™ 1 is defined by Eq. (33). The regularization 



is done in a rotationally invariant way by introducing 
one PV fermion with mass mi and one PV boson with 
mass 

The system of equations for the vertex functions in 
graphical representation can be obtained from that for 
the scalar case, shown in Fig. [7J by changing Sm 2 ,^ —> 
^2,3, #03 -> 903 + Z^mi/j/u-p and by setting 5/j% = 0, 

f3) 

r^J = 0, since we do not consider here the dressing 
of the boson line. After that, the expressions for the 
light-front diagrams are written according to the CLFD 
graph technique rules, taking into account that certain 
lines correspond now to a particle with spin 1/2. 

Expressing through r^ 3 \ as in the scalar case, 
we get a closed matrix equation for the two-body ver- 
tex. Its spin structure is written similarly to Eq. ( [29] ): 

rra/j 



(3) 



u(ki) 



h + b 2 



uj-p 



(58) 



where b\ and b 2 are scalar functions. For clarity, we 
indicate here only the physical components. The PV 
fermion and boson components can be calculated as 
well [15]. 

The three-body component is completely determined 
by four scalar functions, like, for instance 

^(^)r 3 (3) (l,2,3)^(p) (59) 

= *z(fci)(flfi Si + # 2 S 2 + # 3 S 3 + g±S^ju{p), 

where these functions are denoted #1-4, and Si -4 are 
basis spin structures. It is convenient to construct the 
the latter ones as follows 14 : 



Si = 2xi - (1 

m (/) 



uj-p ' 



s 2 = 



uj-p 

iCr> s 



2x\ — (1 — x\ 



m Jj 
1 — — 
uj-p 



75, 



iC, 



m (/) 
uo-p 



-75 



(60) 




Fig. 13 Fermion-boson 3PGF in the three-body approximation. 



with x\ = ^7^, while C ps is the following pseudoscalar: 



m z uj-p 



^k 2fl k 3u p p u 7 . 



(61) 



The function C ps can only be constructed with four 
independent four-vectors. This is the case in LFD for 
n > 3. In the nonrelativistic limit, one would need n > 
4. We can then construct two additional spin structures 
S3 and S4 of the same parity as Si and S 2 by combining 
C ps with parity negative matrices constructed from Si , 
S2, and 75 matrices. 

In our computational procedure we take the limit 
mi — >• 00 analytically and then study the limit /ii —> 00 
numerically. As a result, the calculated vertex functions 
depend parametrically on The main question we are 
interested in concerns the behavior of observables as 
a function of do they remain finite and physically 
reasonable at fii —> 00 or not? For this purpose, we 
calculate the fermion anomalous magnetic moment, us- 



ing the state vector (56). It can be extracted from the 



fermion-boson 3PGF in the three-body approximation, 
as shown in Fig. [l3j The detail of the calculation can 
be found in [15] • We just recall here the main numerical 
results. 

The anomalous magnetic moment is calculated for 
a typical set of physical parameters m = 0.938 GeV, 

H = 0.138 GeV, and two values of the coupling constant 

2 

a = f- =0.2 and 0.5. This mimics, to some extent, a 
physical nucleon coupled to scalar "pions". The typical 
pion-nucleon coupling constant is given by g = 
where (k) is a typical momentum scale, and and 
Fn are the axial coupling constant and the pion decay 
constant, respectively. For (k) — 0.2 GeV we just get 
a -0.2. 

We plot in Fig. [14] the anomalous magnetic moment 

H~~^2 "1 

as a function of log % , for the two different values of 

|_ A* J 

a mentioned above. We show also on each of these plots 
the value of the anomalous magnetic moment calculated 
for the N = 2 truncation, which coincides with the 
anomalous magnetic moment obtained in the second or- 
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Fig. 14 The anomalous magnetic moment in the Yukawa model 
as a function of the PV mass for two different values of the 
coupling constant, a = 0.2 (upper plot) and 0.5 (lower plot). The 
dashed and dotted lines are, respectively, the two- and three- 
body contributions, while the solid line is the total result. The 
value of the anomalous magnetic moment calculated in the N = 2 
approximation is shown by the thin line on the right axis. 



der of perturbation theory. The results for a = 0.2 show 
rather good convergence as \i\ — > oo. The contribution 
of the three-body Fock sector to the anomalous mag- 
netic moment is sizeable but small, indicating that the 



Fock decomposition (16) converges rapidly. This may 



show that once higher Fock components are small, we 
can achieve a practically converging calculation of the 
anomalous magnetic moment. Note that this value of 
a is not particularly small: it is about 30 times the 
electromagnetic coupling, and is about the size of the 
typical pion-nucleon coupling in a nucleus. 

When a increases, we see that the contribution of 
the three-body sector considerably increases. For a = 
0.5 the three-body contribution to the anomalous mag- 
netic moment starts to dominate at large values of fii. 
The dependence of the anomalous magnetic moment 
on the PV boson mass fi\ becomes more appreciable, 
although it keeps rather small. 

In order to have a more physical insight into the 
relative importance of different Fock sectors in the de- 



composition (16) for the state vector, we plot in Fig. 15 



the contributions of the one-, two-, and three-body Fock 
sectors to the norm of the state vector for the two val- 
ues of the coupling constant, considered in this work. 
We see again that at a = 0.2 the three-body contribu- 



0,8 ■ 
0,6 ■ 



0,2 • 
0,0 ■ 



6 7 

log [^ 2 /n 2 ] 



0,8 
0,6 

c 
M 

0,4 
0,2 
0,0 

3456789 10 

log [niVjx 2 ] 

Fig. 15 Individual contributions of the one- (dashed line), two- 
dotted line), and three-body (solid line) Fock sectors to the 
norm, fixed to 1, of the state vector, as a function of the PV 
boson mass /xi, for a = 0.2 (upper plot) and a = 0.5 (lower 
plot). 



tion to the norm is small, while it is not negligible and 
increases with when a = 0.5. 

According to renormalization theory, the PV boson 
mass should be much larger than any intrinsic momen- 
tum scale present in the calculation of physical observ- 
ables. With this limitation, physical observables should 
be independent of any variation of the PV boson mass, 
within an accuracy which can be increased at will. This 
is what we found in our numerical calculation for small 
enough values of a. 

In order to understand the possible origin of the 
residual dependence of the anomalous magnetic mo- 



ment on /ii, we plot in Figs. [16] the two-body spin com- 
ponents bi and 62 calculated at S2 = m 2 , as a func- 
tion of x. As we already mentioned in Sec. |3.3[ the 
spin components for physical particles 61 (#2 = m 2 ) and 
62(^2 = m2 ) should be independent of x in an exact 
calculation. Moreover, 62(^2 = m 2 ) should be zero. It 
is here fixed to zero at a given value of x = x* = — ^— , 
by the adjustment of the constant Z u . We clearly see 
in these plots that the on-shell b\ is not a constant, al- 
though its dependence on x is always weak, while 62 
is not identically zero, although its value is relatively 
smaller than that of b\ for a = 0.2, and starts to be not 
negligible for a = 0.5. A similar situation is observed in 
the scalar case, when we remove "by hands", from the 
state vector, the Fock sector with the antiparticle (see 



the dash-dotted curve in Fig. 10). 
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Fig. 16 The spin components b\ (upper plot) and 62 (lower plot) 
of the two-body vertex function calculated at S2 = m 2 , as a 
function of x, for a = 0.2 (dashed line) and a = 0.5 (solid line), 
for a typical value of [i\ = 100 GeV. 



It is instructive to study the properties of 61,2(52 = 
m 2 ) in perturbation theory. This can be done by calcu- 



lating the amplitudes of the diagrams shown in Fig. 11 
(but for fermions, of course!). Note that in our Fock 
space truncation (|56|), the contribution of the Fbb in- 



termediate state (the left diagram in Fig. 11) is au- 
tomatically taken into account by the solution of the 
eigenstate equation, while the contribution of the FFF 
state (the right diagram in the same figure) is absent. 
If one calculates the sum of both contributions within 
perturbation theory [20], one finds 

m 2 ) = const, 



S2 
S2 



m 2 ) = 0. 



(62) 
(63) 



This is a first indication that the expected properties of 
the on-shell functions b\^ are indeed recovered, when 
antifermion degrees of freedom are involved. This is also 
a confirmation of similar features found in Sec. 14.31 for 
the scalar system. 



masses or electro-weak observables) beyond tree level. 
This is at variance with the meson sector for which a 
meaningful power expansion of any physical amplitude 
can be done. 

While there is not much freedom, thanks to chi- 
ral symmetry, for the construction of the effective La- 
grangian in chiral perturbation theory in terms of the 
pion field — or more precisely in terms of the U field 
defined by U = e %r,7V/ 1 ^ , where f n is the pion decay con- 
stant and r are the Pauli matrices, — one should settle 
an appropriate approximation scheme in order to cal- 
culate baryon properties. Up to now, two main strate- 
gies have been adopted. The first one is to force the 
bare (and hence the physical) nucleon mass to be infi- 
nite, in heavy baryon chiral perturbation theory [21]. In 
this case, by construction, an expansion in characteris- 
tic momenta can be developed. The second one is to use 
a specific regularization scheme [22] in order to separate 
contributions which exhibit a meaningful power expan- 
sion, and hide the other parts in appropriate countert- 
erms. In both cases however, the explicit calculation of 
baryon properties relies on an extra approximation in 
the sense that physical amplitudes are further calcu- 
lated by expanding the effective Lagrangian, denoted 
by C-cpti in a finite number of pion fields. 

Moreover, it has recently been realized that the con- 
tribution of pion- nucleon resonances, like the A and 
Roper resonances, may play an important role in the 
understanding of the nucleon properties at low ener- 
gies [23]. These resonances are just added "by hand" 
in the chiral effective Lagrangian. This is also the case 
for the most important 2ir resonances, like the a and p 
resonances. 

Since in the chiral limit the pion mass is zero, any 
calculation of irN systems demands a relativistic frame- 
work to get, for instance, the right analytical proper- 
ties of the physical amplitudes. The calculation of com- 
pound systems, like a physical nucleon composed of a 
bare nucleon coupled to many pions, relies also on a 
nonperturbative eigenstate equation. While the mass 
of the system can be determined in leading order from 
the iteration of the ttN self-energy calculated in the first 
order of perturbation theory, as indicated in Fig.[l7^a), 
this is in general not possible, in particular, for irN ir- 
reducible contributions, as shown in Fig. [l7^b). 



4.5 Light-front chiral effective field theory 

The calculation of baryon properties within the frame- 
work of chiral perturbation theory is a subject of active 
theoretical developments. Since the nucleon mass is not 
zero in the chiral limit, all momentum scales are a pri- 
ori involved in the calculation of baryon properties (like 




Fig. 17 Iteration of the self-energy contribution in the first order 
of perturbation theory (a); irreducible contribution to the bound 
state equation (b). Dashed lines represent pions. 
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The general framework we have developed above 
is particularly suited to deal with these requirements. 
This leads to the formulation of light-front chiral effec- 
tive field theory (LFxEFT) [24] with a specific effective 
Lagrangian C e f f . The decomposition of the state vector 
in a finite number of Fock components implies to con- 
sider an effective Lagrangian which enables all possi- 
ble elementary couplings between the pion and nucleon 
fields compatible with the Fock space truncation. This 
is indeed easy to achieve in chiral perturbation theory, 
since each derivative of the U field involves one deriva- 
tive of the pion field. In the chiral limit, the chiral effec- 
tive Lagrangian of order p, C P CPT , involves p derivatives 
and therefore at least p degrees of the pion field. In or- 
der to calculate the state vector in the iV-body approx- 
imation, with one fermion and (N — 1) pions, one has 
therefore to include contributions up to 2 (TV — 1) pion 
fields in the effective Lagrangian, as shown in Fig. [THJ 
We thus should calculate the state vector in the TV-body 
approximation with an effective Lagrangian denoted by 
C-^ff and given by 

(64) 



r N _ r p=2(N-l) 
^eff — ^CPT 



(N - 1) pions 



(TV - 1) pions 



Fig. 18 General vertex including a maximum of (N — 1) pion 
fields in the initial and final states. 



While the effective Lagrangian in LF^EFT can be 
mapped out to the CPT Lagrangian of order p, the 
calculation of the state vector does not rely on any mo- 
mentum decomposition. It relies only on an expansion 
in the number of pions in flight at a given light-front 
time. In other words, it relies on an expansion in the 
fluctuation time, ry, of such a contribution. From gen- 
eral arguments, the more particles we have at a given 
light-front time, the smaller the fluctuation time is. At 
low energies, when all processes have characteristic in- 
teraction times larger than tj, this expansion should be 
meaningful. 

It is interesting to illustrate the general features of 
LF%EFT calculations. At order N = 2, we already 
have to deal with irreducible contributions, as shown 



in Fig. 17 ^b). It leads to non-trivial renormalization of 
the coupling constant. The calculation at order N = 3 
incorporates explicitly contributions coming from tttt 
interactions, as well as all low energy 7riV resonances, 
like the A or Roper resonances, as shown schemati- 



the irN state can couple to both J = T = 3/2as well 
as J = T = 1/2 states. We can generate therefore all 
ttN resonances in the intermediate state without the 
need to include them explicitly, provided the effective 
Lagrangian has the right dynamics to generate these 
resonances. This is the case, by the construction, in 
CPT. 



a — p resonances 
A — R resonances 



p(3) 



cally in Fig. (19). Indeed, in the |7T7tA^) Fock sector, 



Fig. 19 Three-body vertex function which exhibits both the con- 
tribution of 27r resonances and A as well as Roper nucleon reso- 
nances. 



Preliminary results obtained by using the PV regu- 
larization scheme can be found in 125 1. 



5 Conclusion and perspectives 

The understanding of the properties of relativistic com- 
pound systems from an elementary Hamiltonian in nu- 
clear and particle physics demands to develop a non- 
perturbative framework. This framework should include 
a well-defined strategy for approximate calculations of 
these properties and a systematic way to improve the 
accuracy. 

We have described in this review a general frame- 
work based on light-front dynamics. In this scheme, the 
state vector of any system of interacting particles is de- 
composed in Fock components. Since for obvious prac- 
tical reasons this decomposition should be truncated to 
take into account a finite number of Fock components, 
we have shown how to control in a systematic way the 
convergence of such expansion. 

Our formalism relies, first, on the covariant formula- 
tion of light-front dynamics, and, second, on a system- 
atic nonperturbative renormalization scheme in order 
to avoid any uncancelled divergences. The applications 
we have presented on QED, on a purely scalar system, 
on the Yukawa model, and on chiral effective field the- 
ory on the light-front have shown the flexibility, the real 
advantages, and the nice features of our formalism. 

Several developments will soon achieve to settle a 
complete framework to deal with field theory on the 
light front. This includes first the full account of an- 
tiparticle degrees of freedom in order to recover, order 
by order in the Fock expansion, the scale invariance 
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of physical observables for arbitrary values of the cou- 
pling constant. This scale invariance should be checked 
within a given regularization scheme. 

We have used up to now the Pauli-Villars regular- 
ization scheme. While this scheme is systematic and 
can be applied to a variety of physical systems, it may 
be cumbersome to implement from a numerical point 
of view in higher order calculations, since it involves 
many non-physical components including Pauli-Villars 
fields. It also demands to perform calculations with very 
large mass scales (the Pauli-Villars fermion and boson 
masses). The use of the Taylor-Lagrange renormaliza- 
tion scheme [26] has proven to be a very natural scheme 
in light-front dynamics. It should be applied to more in- 
volved calculations. 

Finally, it is now well-recognized that the under- 
standing of spontaneous symmetry breaking phenom- 
ena on the light-front can be achieved by taking into ac- 
count nonperturbative zero mode contributions to field 
operators [27]. It remains to include these zero modes 
in the general framework we developed in this review. 
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